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ABSTRACT 

We analyze the distribution of G and K type stars towards the Galactic poles using RAVE and ELODIE radial velocities, 2MASS photometric star 
counts, and UCAC2 proper motions. The combination of photometric and 3D kinematic data allows us to disentangle and describe the vertical 
distribution of dwarfs, sub-giants and giants and their kinematics. 

We identify discontinuities within the kinematics and magnitude counts that separate the thin disk, thick disk and a hotter component. The 
respective scale heights of the thin disk and thick disk are 225± 10 pc and 1048±36 pc. We also constrain the luminosity function and the kinematic 
distribution function. The existence of a kinematic gap between the thin and thick disks is incompatible with the thick disk having formed from 
the thin disk by a continuous process, such as scattering of stars by spiral arms or molecular clouds. Other mechanisms of formation of the thick 
disk such as 'created on the spot' or smoothly 'accreted' remain compatible with our findings. 
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1. Introduction 

It is now widely accepted that the stellar density distribution per- 
pendicular to the Galactic disk traces at least two stellar com- 
ponents, the thin and the thick disks. The change of slope in 
the logarithm of the vertical density distributions at ~ 700 pc 
(Cabr era-Lavers et al. l2005T l or ~ 1500 pc (Gilmore & Reid 
I1983I I above the Galactic plane is usually explained as the signa- 
ture of a transition between these two distinct components: the 
thin and the thick disks. The thick disk is an intermediate stellar 
population between the thin disk and the stellar halo, and was 
initially defined with the other stellar populations by combining 
spatial, kinematic and abundance properties (see a summary of 
the Vatican conference of 1957 by Blaauw ri995l and Gilmore 
& Wyse I1989I I. Its properties are described in a long series of 
publications with often diverging characteristics (see the analy- 
sis by Gilmore [19851 Ojha 2001 Robin et al. 2003 and also by 
Cabrera-Lavers et al. 2005, that give an overview of recent im- 
provements). Majewski ( 1993 ) compared a nearly exhaustive list 
of scenarios that describe many possible formation mechanisms 
for the thick disk. 

In this paper, we attempt to give an answer to the simple 
but still open questions: are the thin and thick disks really two 
distinct components? Is there any continuous transition between 
them? These questions were not fully settled by analysis of star 
counts by Gilmore & Reid ( 1983) and later workers. Other im- 
portant signatures of the thick disk followed from kinematics: 
the age-velocity dispersion relation and also the metallicity- 
velocity dispersion relation. However the identification of a thin- 
thick discontinuity depends on the authors, due to the serious 
difficulty of assigning accurate ages to stars (see Edvardsson et 
al. |1993l and Nordstrom et al. ,2004i ). More recently it was found 
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that the [a/Fe] versus [Fe/H] distribution is related to the kine- 
matics (Fuhrmann |1998t Feltzing et al. 120031 Soubiran & Girard 
120051 Brewer & Carney [20061 Reddy et al. [20061 ) and provides 
an effective way to separate stars from the thin and thick disk 
components. Ages and abundances are important to describe the 
various disk components and to depict the mechanisms of their 
formation. A further complication comes from the recent indi- 
cations of the presence of at least two thick disk components 
with different density distributions, kinematics and abundances 
(Gilmore et al. 2002. Soubiran et al. ,2003; Wyse et al. [2006l) . 

Many of the recent works favor the presently prevailing sce- 
narios of thick disk formation by the accretion of small satellites, 
puffing up the early stellar Galactic disk or tidally disrupting the 
stellar disk (see for example Steinmetz & Navarro [20021 Abadi 
et al. 120031 Brook et al. 120041) . We note however that chemo- 
dynamical models of secular Galactic formation including ex- 
tended ingredients of stellar formation and gas dynamics can 
also explain the formation of a thick disk distinct from the thin 
disk (Samland & Gerhard|2003l Samlandi2004j). 

In this paper, we use the recent RAVE observations of stellar 
radial velocities, combined with star counts and proper motions, 
to recover and model the full 3D distributions of kinematics and 
densities for nearby stellar populations. In a forthcoming study, 
metallicities measured from RAVE observations will be included 
to describe the galactic stellar populations and their history. The 
description of data is given in Sect. 2, the model in Sect. 3, and 
the interpretation and results in Sect. 4. Among these results, 
we identify discontinuities visible both within the density dis- 
tributions and the kinematic distributions. They allow to define 
more precisely the transition between the thin and thick stellar 
Galactic disks. 
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2. Observational data 

Three types of data are used to constrain our Galactic model for 
the stellar kinematics and star counts (the model description is 
given in Se ct. 3): the Two-Micron All-Sky Survey (2MAS S PSC; 
Cutri et al. l2005] l magnitudes, the RAVE (Steinmetz et al. l2006] l 
and ELODIE radial velocities, and the UCAC2 (Zacharias et al. 
I2004I I proper motions. Each sample of stars is selected indepen- 
dently of the other, with its own magnitude limit and coverage 
of sky due to the different source (catalogue) characteristics. 

(1) We select 22 050 2MASS stars within an 8-degree radius 
of the South and North Galactic Poles, with otk magnitudes be- 
tween 5-15.4. Star count histograms for both Galactic poles are 
used to constrain the Galactic model. 

(2) We select 105 170 UCAC2 stars within a radius of 16 
degrees of the Galactic poles, with tok 2MASS magnitudes be- 
tween 6-14. We adjust the model to fit histograms of the yUy and 

proper motion marginal distributions; the histograms com- 
bine stars in 1.0 magnitude intervals for mK=6 to 9 and 0.2 mag- 
nitudes intervals for mK=9 to 14. 

(3) We select 543 RAVE stars ( with 2MASS magni- 
tudes from 8.5 to 1 1.5) within a radius of 15 degrees of the SGP. 
We group them in three histograms according to magnitudes. 
We complete this radial velocity sample with 392 other similar 
stars: TYCHO-II stars selected towards the NGP within an area 
of 720 square degrees, with B-V colors between 0.9-1.1. Their 
magnitudes are brighter than mK=8.5, they were observed with 
the ELODIE spectrograph and were initially used to probe the 
vertical Galactic potential (Bienayme et al. 120061 ). All these ra- 
dial velocity samples play a key role in constraining the vertical 
velocity distributions of stars and the shape of the velocity ellip- 
soid. 
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Fig. 1. Mk/J-K HR Diagram from Hipparcos stars with cr„ln < 
0.1 cross-matched with the 2MASS catalogue. Vertical dashed 
lines represent our color selection J-K= [0.5-0.7] 



2.1. Data selection 

In this paper, we restrict our analysis to stars near the Galactic 
poles with J-K colors between 0.5-0.7 (see figH])- This allows us 
to recover some Galactic properties, avoiding the coupling with 



other Galactic parameters that occurs in other Galactic directions 
(density and kinematic scale lengths, Oort's constants, R^), Vq...). 

The selected J-K=[0.5-0.7] color interval corresponds to 
K3 -K7 d warfs and G3-K1 giants (Koornneef 1983; Ducati et 
al. 1200 11 1. They may be G or K giants within the red clump 
region (the part of the HR diagram populated by high metal- 
licity He-burning core stars). The absolute magnitudes of red 
clump stars are well defined: nearby HIPPARCOS clump stars 
have a mean absolute magnitude Mk - -1.61 with a disper- 
sion of ~ 0.22 (Alves .200 0, see Cannon [19701 for the first 
proposed use of clump stars as distance indicators, see also 
Salaris & Girardi " 20021 Girardi et al. 1998 and other references 
in Cabrera-Lavers et al. 1200 5). This mean absolute magnitude 
does not vary significantly with [Fe/H] in the abundance range 
[-0.6,0] (Alves '2000V Studying nearby stars in 13 open clus- 
ters and 2 globular clusters, Grocholski & Saraiedini 120021 find 
that the mean absolute magnitude of clump stars is not depen- 
dent on metallicity when the [Fe/H] abundance remains within 
the interval [-0.5,0.1]. Saraiedini [20041 finds that, at metalHc- 
ity [Fe/H] =-0.76, the mean absolute magnitude of red clump 
stars drops to Mk - -1.28, a shift of 0.33 mag. Most of the gi- 
ants with metallicity [Fe/H] lower than -0.8 dex are excluded 
by our color selection from our sample. Hence, we did not 
model giants of the metal-weak thick disk, first identified by 
Norris'1985'(see also, Morrison, Flynn & Freeman 1990). This 
represents however only a minor component of the thick disk. 
Although, Chiba & Beers (2000) find that ~ 30 % of the stars 
with -1 > [Fe/H] > -1.7 are thick disk stars, but stars with 
[Fe/H] < -1 represent only 1 per cent of the local thick disk 
stars (Martin & Morrison |19981) . 

K dwarfs within the J-K=[0.5-0.7] color interval also 
have well defined absolute magnitudes that depend slightly 
on metallicity and color We determine their mean ab- 
solute magnitude, Mk=4.15, from nearby HIPPARCOS 
stars using color magnitude data provided by Reid (see 
[http://www-int.stsci.edu/~inr/cmd.html). From Padova 
isochrones (Girardi et al. 120021 ). we find that the absolute 
magnitude varies by 0.4 magnitudes when J-K changes from 
0.5 to 0.7. A change of metallicity of A[Fe/H]=0.6 also changes 
the magnitude by about 0.3, in qualitative agreement with ob- 
served properties of K dwarfs (Reid 1998 , Kotoneva et al 1 20021 ). 
Thus, we estimate that the dispersion of absolute magnitude of 
dwarfs in our Galactic pole sample is ~0. 2-0.4. 

Another important motivation for selecting the J-K=[0.5- 
0.7] color interval is the absolute magnitude step of 6 magni- 
tudes between dwarfs and giants. This separation is the reason 
the magnitude distributions for these two kinds of stars are very 
different towards the Galactic poles. If giants and dwarfs have 
the same density distribution in the disk, in the apparent mag- 
nitude count, giants will appear before and well separated from 
dwarfs. Finally we mention a convenient property of the Galactic 
pole directions: there, the kinematic data are simply related to 
the cardinal velocities relative to the local standard of rest (LSR). 
UCAC2 proper motions are nearly parallel to the U and V ve- 
locities, and RAVE radial velocities are close to the vertical W 
velocity component. 

2.2. How accurate is tlie available data? 

The star magnitudes are taken from the 2MASS survey which is 
presently the most accurate photometric all sky survey for prob- 
ing the Galactic stellar populations. Nevertheless, since our color 
rang is narrow, we have to take care that the photometric errors 
on J and K do not bias our analysis. 
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The mean photometric accuracy ranges from 0.02 in K and 
J at magnitudes mK=5.0, to 0.15 in K and 0.08 in J at magni- 
tudes mK=15.4. The error in J-K is not small considering the 
size (A(J-K) = 0.2) of the analyzed J-K interval, 0.5 to 0.7. We 
do not expect, however, that it substantially biases our analy- 
sis. For wk brighter than 10, the peak of giants is clearly iden- 
tified in the J-K distribution within the J-K= [0.5-0.7] interval 
(see Fig|2]or Figure 6 from Cabrera-Lavers et al., 2005). This 
peak vanishes only beyond mj^-ll. At fainter K magnitudes, the 
dwarfs dominate and the J-K histogram of colors has a constant 
slope. This implies that the eiTor in color at faint magnitudes 
does not affect to first order the star counts. 

We find from the shape of the count histograms that, in the 
direction of the Galactic Pole and with our color selection J-K - 
[0.5-0.7] the limit of completeness is ~15.5-15.6. Moreover, 
the contamination by galaxies must be low within the 2MASS 
PSC. It is also unlikely that compact or unresolved galaxies are 
present: according to recent deep J and K photometric counts 
(see Figure 15 of lovino et al. 120051 ), with our color selection, 
galaxies contribute only beyond otk ~16. We conclude that we 
have a complete sample of stars for magnitudes from 5.0 to 15.4 
in K, towards the Galactic poles. 

The UCAC2 and RAVE catalogues however are not com- 
plete. Making it necessary to scale the proper motions and ra- 
dial velocities distributions predicted by our model for complete 
samples. The total number of stars given by the model for the 
distribution of proper motions (or radial velocity) in a magni- 
tude interval is multiplied by the ratio between the number of 
stars observed in UCAC2 (or RAVE) divided by the number of 
stars observed in 2MASS. 
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Fig. 2. K / J-K Color Magnitude Diagram obtained with 
2MASS stars within a 8 degrees radius around the North 
Galactic pole. Dashed lines represent the limit of our color se- 
lection: J - K = [0.5,0.7]. 

Towards the North Galactic Pole (NGP), the error on 
the UCAC2 proper motions used in our analysis varies from 
Imasyr"' for the brightest stars to 6masyr"' at mK=14. 
Towards the South Galactic Pole (SGP), the error distribution 
looks similar, with the exception of a small fraction of stars with 
mK from 11 to 14 having errors around 8 or 13 masyr"'. The 



only noticeable difference between the histograms at the NGP 
and SGP is that the peak of the proper motion distribution is 
slightly more flattened at the SGP, for magnitudes >13 (see 
fig- El). This difference is related to the different eiTor distribu- 
tions towards the NGP and SGP. 

The analyzed stars are located at distances from 200 pc to 
1 kpc for dwarfs and to 1 .5 kpc for giants. A 2 mas yr"' error rep- 
resents lOkms"' at 1 kpc, and 6 mas yr"', an error of 30km s"'. 
This can be compared to the era values for the isothermal com- 
ponents, for instance ~60kms ' for the thick disk that is the 
dominant stellar population 1.5 kpc from the plane. Adding the 
errors in quadrature to the velocity dispersion would modify a 
real proper motion dispersion of 60kms ' to an apparent dis- 
persion of 67kms ^ The apparent dispersion would be only 
60.8 km s"' if the stars have a 2masyr"' accuracy. Therefore, 
we overestimate the c^^ dispersion of the thick disk by 5 to 10 
percent. This eff'ect is lower for the thin disk components (the 
stars are closer and their apparent proper motion distributions 
are broader). We have not yet included the eff'ect of proper mo- 
tion errors within our model. This error has just an impact of the 
determination on the velocity dispersions era and o-y and on the 
ellipsoid axis ratio ctu/ctw of each stellar disk component, but 
does not change the determination of vertical velocity disper- 
sions cTw which are mainly constrained by the magnitude star 
count and the radial velocities. Hence, it is not significant in our 
kinematic decomposition of the Galactic disk. 

The accuracy of proper motions can also be gauged from the 
stability of the peaks of proper motion distributions: comparing 
1 12 /ij/ and fiy histograms for different magnitude intervals, we 
find no fluctuations larger than 3-5 mas yr ' . 

A more complete test is performed by comparing the 
UCAC2 proper motions (with our J-K color selection) to the 
recent PM2000 catalogue (Ducourant et al. 1200 6!) in an area of 
8x16 degrees around ci;2{)(K)=12h50m, iSaooo - 14deg close to 
the NGP. PM2000 proper motions are more accurate, with er- 
rors from 1 to 4 mas yr The mean differences between proper 
motions from both catalogues versus magnitudes and equato- 
rial coordinates do not show significant shifts, just fluctuations 
of the order of ~0.2 masyr We also find that the disper- 
sions of proper motion differences are ~2 mas yr"' for otk <10, 
4 masyr"' with mK=10-13, and 6 masyr"' with mK=13-14. 
These dispersions are dominated by the UCAC2 errors. 

From the internal and external error analysis, RAVE radial 
velocities show a mean accuracy of 2.3 km s"' (Steinmetz et 
al. [2006 ). Radial velocities of stars observed with the ELODIE 
echelle spectrograph are an order of magnitude more accurate. 
These eiTors have no impact on the determination of the vertical 
velocity dispersion of stellar components that ranges from 10 to 
50 km s"' , but the reduced size of our radial velocity samples to- 
wards the poles (about 1000 stars) limits the accuracy achieved 
in modeling the vertical velocity dispersions. 

3. Model of the stellar Galactic disks 

The basic ingredients of our Galactic model are taken from tradi- 
tional works on star count and kinematic modeling, for instance 
see Pritchet (1983 ); Bahcall (1984); Robin & Creze (1986 ). It is 
also similar to the recent developments by Girardi et al. (120051 ) 
or by Vallenari et al. (120061 ). 

The kinematic modeling is entirely taken from Ratnatunga 
et al. (1989) and is also similar to Gould's (2003) analysis. Both 
propose closed-form expressions for velocity projections; the 
dynamical consistency is similar to Bienayme et al. d 1 9871 ) and 
Robin et al. (120031 120041 ). 
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Our analysis, limited to the Galactic poles, is based on a set 
of 20 stellar disk components. The distribution function of each 
component or stellar disk is built from three elementary func- 
tions describing the vertical density p,- (dynamically self con- 
sistent with the vertical gravitational potential), the kinematic 
distribution (3D-gaussians) and the luminosity function 0,^. 

We define N{z, Vr,V^,V,; M) to be the density of stars in the 
Galactic position-velocity-(absolute magnitude) space 

^ = A<z)./;(V«, V^, Vc)0,7t(M) 
ik 

the index / differentiates the stellar disk components and the 
index k the absolute magnitudes used to model the luminosity 
function. 

From this model, we apply the generalized equation of stellar 
statistics: 

A(m, ^iu yu;„ K) = J N(z, V«, V^, K; M) oj dz 

to determine the Aim) apparent magnitude star count equation 
as well as the marginal distributions of both components /i; and 
lib of proper motions and the distributions of radial velocities for 
any direction and apparent magnitudes. For the Galactic poles, 
we define fiu and fiy as the proper motion components parallel to 
the cardinal directions of U and V velocities. For a more general 
inverse method of the equation of the stellar statistic, see Pichon 
et al. ( [2002] i. 

3.1. The vertical density 

Each stellar disk is modeled with an isothermal velocity distribu- 
tion, assuming that the vertical density distribution (normalized 
at z=0) is given by the relation: 

p,(z)=exp(-0(z)/4,) (1) 

where <l'(z) is the vertical gravitational potential at the solar 
Galactic position and cr.j , is the vertical velocity dispersion of 
the considered stellar component /. The Sun's position Ze above 
the Galactic plane is also used as a model parameter Such ex- 
pressions were introduced by Oort d 19221 1, assuming the station- 
arity of the density distributions. They ensure the consistency 
between the vertical velocity and density distributions. For the 
vertical gravitational potential we use the recent determination 
obtained by Bienayme et al. ( 12006! ) based on the analysis of 
HIPPARCOS and TYCHO-II red clump giants. The vertical po- 
tential is defined at the solar position by: 

<D(z) = 4nG (lo ( VJTd2 - d) + peff z') 

with Eo = 48 Mopc"2, D = 800 pc and petf = 0.07 Mgpc^l 

It is quite similar to the potential determined by Kuijken & 
Gilmore (.1989 J and Hohnberg & Flynn (.2004 J . 

3.2. The kinematic distributions 

The kinematical model is given by shifted 3D gaussian veloc- 
ity ellipsoids. The three components of mean streaming motion 
{(U), (V), {W}) and velocity dispersions (ctrr, cr^^, cr^^, referred 
to the cardinal directions of the Galactic coordinate frame, pro- 
vide a set of six kinematic quantities. The mean stream motion is 
relative to the LSR. The Sun's velocity Uq and are model pa- 
rameters. We define the {V) stream motion as: (V) - - V© - Viag. 



We adopt an asymmetric drift proportional to the square of ctrr: 
^lag = ct^rIK, where the coefficient kg is also a model parame- 
ter We assume null stream motions for the other velocity com- 
ponents, thus {U) - -Uq and (W) - -Wq. 

For simplicity, we have assumed that the ctrr/o-^^ ratio is the 
same for all the components. It is well known that the assump- 
tions of a constant ctrr/o-^^ ratio, of a linear asymmetric drift 
and of 2D gaussian U and V velocity distributions hold only for 
cold stellar populations (see for instance Bienayme & Sechaud 
11997 ). These simple assumptions allow a direct comparison with 
similar studies. It allows also an exact integration of count equa- 
tions along the line of sight. Thus the convergence of parameters 
for any single model is achieved in a reasonable amount of time 
(one week). The model includes 20 isothermal components with 
cr,7 from 3.5 to 70kms"'. We choose a step of 3.5 kms"' which 
is sufficient to give a realistic kinematic decomposition and per- 
mit calculation in a reasonable time. The two first components 
(T^z- 3.5 and 7kms"' were suppressed since they do not con- 
tribute significantly to counts for mK>6 and are not constrained 
by our adjustments. The components between 10 and 60 km s"' 
are constrained by star counts, proper motions histograms up to 
magnitude 14 in K and radial velocity histograms for magni- 
tudes mK=[5.5-11.5]). The model includes isothermal compo- 
nents from 60 to 70kms"' to properly fit the star counts at the 
faintest apparent magnitudes > 15.0. All the values of the 
kinematic components depend on the adopted galactic potential. 

The velocity ellipsoids are inclined along the Galactic merid- 
ian plane. The main axis of velocity ellipsoids are set parallel to 
con-focal hyperboloids as in Stackel potentials. We set the focus 
at z/,yp=6kpc on the main axis giving them realistic orientations 
(see Bienayme 1999). The non-zero inclination implies that the 
vertical density distributions of each isothermal component is 
not fully dynamically consistent with the potential. Since the z- 
distances are below 1.5 kpc for the majority of stars with kine- 
matic data, and since the main topic of this paper is not the de- 
termination of the Galactic potential, we do not develop a more 
consistent dynamical model. 

3.3. The luminosity functions 

The luminosity function of each stellar disk component is mod- 
eled with n different kinds of stars according to their absolute 
magnitude: 

(PiiM) = 2_j ^»(^) = -1= Zj '^ike''^~> 

k=\,n V27rO-M A.= i,„ 

where Cik is the density for each type of star (index k) of each 
stellar disk component (index /). 

We use four types of stars to model the local luminosity 
function (see Fig. [3]). More details on the way that we have 
determined it is given in section l4~4l Stars with a mean abso- 
lute magnitude Mk = -1.61 are identified to be the red clump 
giants {k - \) that we will call 'giants', with Mk = -0.89 
and Mk - -0.17 for first ascent giants that we categorize as 
'sub-giants' {k - 2 — J) and Mk = 4.15 are labelled dwarfs 
{k - 4) (see fig. [T]i. We neglected 'sub-giant' populations hav- 
ing absolute magnitude Mk between 0.2 and 2. Their presences 
marginally change the ratio of giants to dwarfs, since their mag- 
nitudes are lower, and their total number in the magnitude counts 
appears significantly smaller than the other components. In fact, 
we initially tried to introduce 10 types of stars (spaced by 0.7 
absolute magnitude intervals). This still improves the fit to the 
data. However due to the small contribution of the 'sub-giants' 
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components with Mk — [0.2 - 2], they were not determined with 
a useful accuracy. We adopt ctm - 0.25, justified by the narrow 
range of absolute magnitudes both for red clump giants and for 
dwarfs on the luminosity function. 

The 4x20 coefficients c,; are parameters of the model. In or- 
der to obtain a realistic luminosity function, we have added con- 
straints to the minimization procedure. For each kinematic com- 
ponent /, we impose conditions on the proportion of dwarfs, gi- 
ants and sub-giants following the local luminosity function. We 
have modeled our determination of the local luminosity function 
of nearby stars (see Fig.[3]l. We obtained : 

- a ratio of the density of dwarfs (k=4) to the density of giants 
(^=1) of 12.0, so we impose: ^ > 10 

- a ratio of the density of giants (^=1) to the density of sub- 
giants(^=2) of 2.3, so we impose: > 2 

- and the density of sub-giants (^=2) is greater than the density 
of sub-giants (A:=3), so we impose: c,j > Ci 2- 

If we do not include these constraints, the various components 
are populated either only with dwarfs or only with giants. 
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Fig. 3. Local luminosity function: The histogram is our deter- 
mination of the local luminosity function for nearby stars with 
error bars. The red (or dark grey) dashed line is a fit of the lu- 
minosity function with four gaussians (blue or light grey line) 
corresponding to the dwarfs, the giants and the two types of sub- 
giants. 



4. Results and discussion 

The 181 free model parameters are adjusted through simulations. 
Each simulation is compared to histograms of counts, proper 
motions and radial velocities (see Sect. |2]for the description of 
data histograms and see Figures H] |5] |6] [T] E]) for the compari- 
son of the best fit model with data. The adjustment is done by 
minimizing a function using the MINUIT software (James 
[2004). Equal weight is given to each of the four types of data 
(magnitude counts, /ijy proper motions, jUy proper motions, and 
radial velocities). This gives relatively more weight to the ra- 
dial velocity data whose contribution in number is two orders of 
magnitude smaller than for the photometry and proper motions. 

By adjusting our Galactic model, we derive the respective 
contributions of dwarfs and giants, and of thin and thick disks. 
One noticeable result is the kinematic gap between the thin and 



thick disk components of our Galaxy. This discontinuity must be 
the consequence of some specific process of formation for these 
Galactic components. 

Fitting a multi -parameter model to a large data-set raises the 
question of the uniqueness of the best fit model, and the robust- 
ness of our solution and conclusions. For this purpose, we have 
explored the strength of the best Galactic model, by fitting vari- 
ous subsets of data, by modifying various model parameters and 
adjusting the others. This is a simple, but we expect efficient, 
way to understand the impact of parameter correlations and to 
see what is really constrained by model or by data. A summary 
of the main outcomes is given below. 

From these explorations, we choose to fix or bound some 
important Galactic model parameters which would otherwise be 
poorly constrained: i) we fix the vertical Galactic potential (ad- 
justing the K- force does not give more accurate results than for 
instance in Bienayme et al., 2006 since we only increase by a 
factor 2 the number of stars with measured radial velocities), 
ii) the asymmetric drifts of all kinematic components are linked 
through a unique linear asymmetric drift relation with just one 
free parameter; the solar velocity component Vq is also fixed, iii) 
the axis ratio of the velocity ellipsoids is bounded; for thin disk 
components (crw < 25kms"') we set cru/crw > 1.5, for thick 
disks (cTw > 30kms"', o-u/ctw > 1.1). 

The agreement between our fitted model and the observed 
counts is illustrated by the various magnitude, proper motion and 
radial velocity distributions (Figs.|4]|5]|6]|2l[8j. We can consider 
that globally the agreement is good, if we note the small val- 
ues obtained. We just comment the main disagreements visible 
within these distributions. They can be compared to recent sim- 
ilar studies (Girardi et al. 120051 Vallenari et al. 2006). 

The agreement for the apparent magnitude distribution looks 
satisfying in Fig.|4] 

The comparison of observed and modeled fiu proper motion 
distributions does not show satisfactory agreement close to the 
maxima of histograms at apparent magnitude otk <10 ( NGP or 
SGP, see Fig. |5]l. We have not been able to determine if this is 
due to the inability of our model to describe the observed data, 
for instance due to simplifying assumptions (gaussianity of the 
velocity distribution, asymmetric drift relation, constant ratio of 
velocity dispersions, etc.). We note that this disagreement may 
just result from an underestimate of the impact of the proper 
motion errors. 

Some possible substructures are seen in proper motion his- 
tograms for the brightest bins (otk <7, Fig. |5]i; they are close 
to the level of Poissonian fluctuations and marginally signifi- 
cant. One of the possible structures corresponds to the known 
Hercules stream (U = -42 kms"' and V - -52kms"', Famaey 
et al 2005 ). 

For faint magnitude (wk >1 1) bins (Fig.|6]l, small shifts (~ 3- 
5 mas yr"') of fiu explain most of the diff'erences between North 
and South and the larger 

At hjk within 10-13 (Fig. |6]l, the wings of jiu histograms 
look slightly different between North and South directions; it 
apparently results from shifts of North histograms versus South 
ones. 

A disagreement of the model versus observations also ap- 
pears within the wings of fiy distributions, (wk within 10-13, 
Fig. |6]l. This may introduce some doubt concerning our ability 
to correctly recover the asymmetric drift, because the negative 
proper motion tail of fiy distributions directly reflects the asym- 
metric drift of the V velocity component. However, we estimate 
that our determination of the asymmetric drift coefficient is ro- 
bust and marginally correlated to the other model parameters. 
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Fig. 4. Magnitude count histogram towards the North Galactic Pole. Left: model prediction (dashed line) is split according to star 
types: giants (red or black line), sub-giants (dot-dashed and dotted) and dwarfs (green or grey line). The right figure highlights the 
contributions of thin and thick disks (respectively thin and thick lines), for dwarfs (green or grey) and giants (red or black). 
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These comparisons of observed and model distributions sug- 
gest new directions to analyze data. In the future, we plan to use 
the present galactic model to simultaneously fit the RAVE ra- 
dial velocity distribution in all available galactic directions. This 
result will be compared to a fit of our model to proper motion 
distributions over all galactic directions. This will give a better 
insight into the inconsistency between radial velocity and proper 



motion data, and also for possible inconsistency in our galactic 
modeling. 

4. 7 . The transition from dwarfs to giants 

Within the J-K=[0.5-0.7] interval, the proper motion is an excel- 
lent distance indicator: there is a factor of 14 between the proper 
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motion of a dwarf and the proper motion of a giant with the same 
apparent magnitudes and velocities. Combining proper motions 
and apparent magnitudes, our best-fit Galactic model allows us 
to separate the contributions of dwarfs and giants (Fig.|4|i. 

We deduce that, towards the Galactic poles, most of the 
bright stars are giants. At = 7.2, only 10% are dwarfs and 
at wk = 9.6 only 50% are giants. We have checked if the con- 
tribution of sub-giants with absolute magnitude Mk - [0.2 - 2] 
can change the contribution of dwarfs and giants. At < 10, 
the contribution of sub-gaint with Mk = [0.2 - 2] is at least one 
order of magnitude lower So the ratio of giants and dwarfs is 
unchanged. Furthermore, the RAVE data confirm our model pre- 
diction. This is in contradiction with Cabrera-Lavers et al. (I2005I I 
statement based on the Wainscoat et al. (|1992|) model which esti- 
mates that, at magnitude < 10, giants represent more than 90 
% of the stars. The Wainscoat model assumes only one disk with 
a scale height of 270 pc for the giants and 325 pc for the dwarfs. 
In our model, we find a scale height of 225 pc both for the giants 
and the dwarfs. This explains why we find more dwarfs at bright 
magnitudes ( < 10). 

Faint stars are mainly dwarfs, 80% at = 11.6 while at 
otk - 11.9, only 10% are giants. The 50%-50% transition be- 



tween giants-sub-giants and dwarfs occurs at mK ~ 10.1. This is 
a robust result from our study that depends slightly on the abso- 
lute magnitude adopted for dwarf and giant stars. We have not 
tried to change our color range. If we take a broader color inter- 
val, the dispersion around the absolute magnitude of dwarfs will 
be larger, but our results are not expected to change. For another 
color interval, we can expect this result to be different, since we 
would be looking at a different spectral type of star. 
A confirmation of the dwarf-giant separation between magni- 
tudes otk = [5.5 - 11.5] comes from RAVE spectra. With the 
preliminary determination of the stellar parameters (Tpyy, \og{g) 
and [Fe/H]) of RAVE stars, we choose to define giant stars with 
log(,g) < 3 and dwarfs with \og(g) > 4. The comparison of the 
number of giants and dwarfs predicted by our best model to the 
observed one is in good agreement (see fig. [8]). 

4.2. The scale heights of stellar components 

Our dynamical modeling of star counts allows us to recover the 
vertical density distribution of each kinematic component p,(z), 
with the exact shapes depending on the adopted vertical poten- 
tial <l>(z). We recover the well-known double-exponential shape 
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Distance (pc) 

Fig. 9. Model of the vertical stellar density p{z) towards the 
the North Galactic Pole (dashed line) and its thin and thick 
disk decomposition (respectively thin and thick lines). The 
thin disk includes the isothermal kinematic components with 
(Tw <25kms the thick disks include components with crw > 
25kms-i. 



of the total vertical number density distribution ptot(z) (Fig. |9l). 
Since we estimate that the kinematic decomposition in isother- 
mal components is closer to the idealized concept of stellar pop- 
ulations and disks, we identify the thin disk as the components 
with vertical velocity dispersions o-^ smaller than 25kms"' 
and the thick disk with crw from 30 to 45.5 km s"' (Fig. [T2b . 
Following this identification, we can fit an exponential on the 
thin and thick disk vertical density component (thin line and 
thick lines respectively of Fig. |9]l. The scale height of the thin 
disk is 225±10pc within 200-800pc. For the thick disk, within 
0.2-1.5 kpc, the scale height is 1048±36pc. If we consider all the 
kinematic components without distinguishing between the thin 
and thick disk, we can fit a double exponential with a scale length 
of the thin disk217±15pc and of the thick disk 1064+38pc. We 
calculate the error of the scale length from the error on the in- 
dividual kinematic disk components 0^,,, / (see Tab. [T]i. We have 
performed a Monte-Carlo simulation on the value of the compo- 
nents and obtained the error bars for the scale length of the thin 
and thick disk both independently and together. 

We note that our density distribution is not exponential for 
z < 200 pc: this mainly results from the fact that we do not model 
components with small velocity dispersions crw <8kms 
Thus our estimated density at z=0 cannot be directly compared, 
for instance, to Cabrera-Lavers et al. ( I2005I I results. With this 
proviso, the star number density ratio of thick to thin disk stars 
at z=Opc is 8.7% for the dwarfs. 

One candidate to trace the thin and thick disk are the red 
clump giants. In fact, at z-distances larger than ~500pc (i.e. 
otk larger than ~7.0, see Fig. 21 there are more thick disk gi- 
ants than thin disk giants. Cabrera-Lavers et al. ( 120051 ) have an- 
alyzed them using 2MASS data. To do this, they select all stars 
with color J-K= [0.5-0.7] and magnitude < 10. But, beyond 
magnitude 9, the proportion of giants relative to sub-giants and 
dwarfs decreases quickly. At mi^-9.6, giants just represent half 
of the stars, and their distance is about 1 .7 kpc. Thus, we must be 
cautious when probing the thick disk with clump giants and we 
have first to determine the respective sub-giant and dwarf contri- 
butions. However, Cabrera-Lavers et al. ([2005) obtained a scale 
height of 267±13pc and 1062+52pc for the thin and thick disks 
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which is in relatively good agreement with the values obtained 
from our model. 

For dwarfs that dominate the counts at faint apparent mag- 
nitudes otk > II (distances larger than ~ 240 pc), we use the 
photometric distance: 

ZpHo, = 10^""^-^'^-'^' (2) 

where Mk is equal to 4.15 (the value for the dwarfs). 

Doing so, we obtain the number density n{zph„i) of stars seen 
along the line of sight at the SGP and NOP (Fig. [TO]). These 
plots show a well-defined first maximum at Zp/,o,=500pc (SGP) 
or 700 pc (NGP) related to the distribution of thin disk dwarfs. 
At 0.9-1.1 kpc, n{zphot) has a minimum and then rises again at 
larger distances, indicating the thick disk dwarf contribution. 

However, the use of photometric distances can introduce a 
systematic error for thick disk dwarfs that have lower metal- 
licities . The mean metallicity of the thick disk population at I 
kpc is <[Fe/H]) ^ -0.6 (Gilmore et al. [19951 Carraro et al. [19981 
Soubiran et al. 2003). 

The metallicity variation from [Fe/H]=0.0 for the thin disk to 
[Fe/H]=-0.6 for the thick disk means that the absolute magnitude 
Mk changes from 4. 15 to 4.5. So, we smoothly vary the absolute 
magnitude with the metallicity from the thin to the thick disk, in 
this way: 

MK([f e///]) = Mk,o + 0.035mK (3) 

where Mk,o is equal to 4. 15. 

The counts continue to show two maxima (Fig.[TTI). even if 
the minimum is less deep. The minimum delineates a discontin- 
uous transition between the thin and thick components. 

The superposition of the model on the number density 
n(Zphoi) shows only approximate agreement (Fig.lTOb. We think 
that is due to non-isothermality of the real stellar components. 
Anyway, the fact that the model does not reproduce exactly the 
observation does not weaken the conclusion about the kinematic 
separation of the thin and thick disk. It reinforces the need 
for a clear kinematic separation between the two disks in the 
kinematic decomposition (Fig.[T2ll. 

We also notice, in Fig. [TOl the difference in counts between 
the North and the South. This difference allows us to deter- 
mine the distance of the Sun above the Galactic plane, Zq = 
-1-20.0 + 2.0 pc, assuming symmetry between North and South. 
We also note that the transition between thin and thick disks is 
more visible towards the SGP than towards the NGP. 

4.3. The thin-thick disk transition, and the kinematic 
distribution function 

The minimum at z ~ 1 kpc in the n{z) distribution (Fig. [TOl i 
provides very direct evidence of the discontinuity between 
stellar components with small velocity dispersions {crw-lO- 
25kms"') and those with intermediate velocity dispersions 
(crw -45.5 kms-i) (left panel Fig.inii. 

Another manifestation of this transition is well known from 
the logp(z) density distribution (Fig.|9l) which shows a change of 
slope at z=500-700pc. This feature can be successfully modeled 
with two (thin and thick) components (e.g. Reid and Gilmore 
119831 ), which is an indication of a discontinuity between the thin 
and thick disks of our Galaxy. 

It is conclusive evidence, only if we show that we can not fit 
accurately the star counts or vertical density distributions with a 
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Fig. 10. Data (histogram with error bars) and model (dashed hne) for the NGP (left) and SGP (right) vertical density distribution 
using photometric distances nphoti^) for dwarf stars. The transition between thin and thick components is revealed by a minimum at 
z ~1 kpc. The main contributing components are plotted, for the thin disk (thin continuous line) ctw = 10.5 (dot-dashed), 14 & 17.5 
(triple dot-dashed), 21 & 24.5 km s~' (dotted) and for the thick disk (thick continuous line) (tw = 45.5 km s ' . 




Fig. 11. Histograms of the vertical density distribution for the NGP (left) and SGP (right) using photometric distances riphoti^) for 
dwarf stars with a smooth variation in the [Fe/H] from the thin to the thick disk. 




Fig. 12. Left: The local ctw kinematic distribution function. The contributing components to star counts can be put together in a thin 

disk component {cr^/^r <25 km s~'), a thick disk (isothermal with cry^-A5.5 km s"^) and a hotter component with ctw ~65 km s"' . The 
two first components with crw-^-i and 7 kms"^ are set to zero by construction. Right: A Kinematic Distribution Function (KDF) 
that tries to reproduce the magnitude star counts and the kinematic data: this model has been obtained requiring the continuity of 
the KDF from cr„= 10 to 48 km s'^ 
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Table 1. List of the values of the kinematic disk components 
(pkinj ( 10''x number of stars / pc-' ) with the individual errors 
absolutes and relatives in percent. 



continuous set of kinematic components (without a gap between 
the thin and the thick disks). We find that the constraint of a set of 
kinematic components following a continuous trend (right panel 
of Fig.fTSli raises the reduced;^'^, in particular on SGP magnitude 
counts, from 1.59 to 3.40. This confirms the robustness of our 
result and conclusion on the wide transition between thin and 
thick stellar disk components. 

Adjusting the Galactic model to star counts, tangential and 
radial velocities, we can recover the details of the kinematics 
of stellar populations, and we determine the local ctw kinematic 
distribution function (left panel of Fig.[T2land Tab.[T]). This kine- 
matic distribution function clearly shows a large step between 
the kinematic properties of the thin and thick disks. We define 
the thin disk as the components with crw covering 10-25 km s"', 
and the thick disk as the components with ctw covering 30-45 
kms"'. The counts and radial velocities by themselves already 
show the kinematic transition that we obtain in the kinematic de- 
composition. The fit of proper motions confirms the conclusion 
from the star counts and radial velocities, even if a fraction of the 
proper motions jU/ and fit> at magnitude mK fainter than 13 have 
significant errors (> 20 kms '). The only consequence for the 
proper motion errors is that we obtained an ellipsoid axis ratio 
o'u/o'w different from the classical values (see Sec. 14.51 ). 

The last non-null components at approximately crw ~ 
65 km s"' are necessary to fit the faintest star counts at ~ 15. 
But, they do not result from the fit of proper motion histograms 
(since, unfortunately, they stop at ~14). Thus their exact na- 
ture, a second thick disk or halo (they would have very different 
asymmetric drift) cannot be solved in the context of our analysis. 

4.4. The luminosity function of stellar components 

Our distant star count and kinematic adjustment constrains the 
local luminosity function (LF). We make the comparison with 
the local LF determined with nearby stars. However, the bright- 
est HIPPARCOS stars needed to determine the local LF are satu- 
rated within 2MASS and have less accurate photometry. We can 



also compare it to the LF determined by Cabrera-Lavers et al. 
(12005 ") who use a cross-match of HIPPARCOS and MSX stars 
and estimate magnitudes from MSX A band magnitudes 
(hereafter [8.3]). However we note from our own cross-match 
of HIPPARCOS -MSX-2M ASS (non saturated) stars that their 
LF, for stars selected from V-[8.3], corresponds mainly to stars 
with J-K colors between 0.6-0.7 rather than between 0.5-0.7. A 
second limitation for a comparison of LFs is that our modeling 
does not include the stellar populations with small velocity dis- 
persions (cTjy < 8 kms"'). For these reasons, we determine a 
rough local LF based on 2MASS-HIPPARC0S cross-matches, 
keeping stars with V<7.3 or distances <125pc, and using the 
color selection V-K between 2.0 and 2.6, that corresponds ap- 
proximately to J-K - [0.5-0.7]. Using V and K magnitudes min- 
imizes the effects of the J-K uncertainties. Considering these 
limitations, there is reasonable agreement between the local LF 
obtained with our model using distant stars and the LF obtained 
from nearby Hipparcos stars (see Fig.[T3]l. 
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Mk 

Fig. 13. The local Luminosity Function of K stars from our mod- 
eling of star counts towards the Galactic poles (line) compared 
to the LF function from nearby Hipparcos K stars by Cabrera- 
Lavers et al. (2005) (red or black histogram) and our own es- 
timate of the local LF: see text (green or grey histogram with 
error bars). The scale of Cabrera et al.'s LF has been arbitrarily 
shifted. 



4.5. The Stellar kinematics 

Many of the stellar disk kinematic properties obtained with our 
best fit Galactic model are comparable with previously pub- 
lished results. We make the comparison with the analysis of 
HIPPARCOS da ta (De hnen & Binney [T9981 Bienayme [19991 
Nordstrom et al. 120041 Cubarsi & Alcobe 2004; Famaey et al. 
2005), and also with results published from remote stellar sam- 
ples using a wide variety of processes to identify thin and thick 
kinematic components (Bartasiute 1994, Flynn & Morrel [79971 
Soubiran et al. 2003; Pauli et al. 2005). 

We obtain for the Sun motion relative to the LSR, 
Mo=8. 5+0.3 kms"' and Wo^ll.lil.Okms"'. We find for the 
asymmetric drift coefficient, fcfl=76±4kms"', compared to 
80±5kms"' for nearby HIPPARCOS stars (Dehnen & Binney 
1998) and the thick disk lag is Vi^g = o-^/Ar^ = 33±2kms"' rel- 
ative to the LSR. We note that this value of the thick disk lag is 
close to the value of Chiba & Beers (,2000) and other estimates 
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prior to this. It is less in agreement with the often-mentioned val- 
ues of 50-lOOkms ' from pencil-beam samples. These may be 
more affected by Arcturus group stars which are more dominant 
at higher z-values. 

Our determination of the asymmetric drift coefficient is 
highly correlated to V©. The reason is that we do not fit pop- 
ulations with low velocity dispersions and small Viag since we 
do not fit star counts with < 6: as a consequence the slope 
of the relation, Viag versus cru, is less well constrained. To im- 
prove the ka determination, we adopt = 5.2 km s"' (Dehnen 
& Binnev ll998t Bienayme 19991) . The adjusted cru/o-y veloc- 
ity dispersion ratio, taken to be the same for all components, is 
1.44 ± 0.02. We obtain cru/crw ratios significantly smaller than 
those published using nearby samples of stars. For the thin disk 
components, we find cr(//crjy=1.50 to 1.62 (compared to pub- 
lished values ~2 by authors using HIPPARCOS stars). For the 
thick disk, we obtain crf//crjy=l.l, instead of ~ 1.5-1.7 typi- 
cally obtained with nearby thick disk stars by other authors. 

While there is no dynamical reason preventing the varia- 
tion of (Tu/o-w with z, we suspect that our low cru/crw ratio 
at large z for the thick disk results from a bias within our model 
due to the outer part of the wings of some proper motion his- 
tograms not being accurately adjusted. This may be the conse- 
quence of an incorrect adopted vertical potential or, as we think, 
more likely the non-isothermality of the real velocity distribu- 
tions. This suspicion is reinforced since fitting each proper mo- 
tion histograms separately with a set of gaussians gives us larger 
values for cra/crw- 

Our results can be directly compared with the very recent 
analysis by Vallenari et al. (2006) of stellar populations towards 
the NGP using BVR photometry and proper motions (Spagna 
et al. 119961 ). Their model is dynamically consistent but based 
on quite different hypotheses from ours; for each stellar popula- 
tion, they assume that in the Galactic plane cr^, is proportional to 
the stellar density p (Kruit & Searle 119821 ). They also assume 
that both velocity dispersions, cr|, and crj^g, follow exponen- 
tial laws with the same scale exponential profile as the surface 
mass density (Lewis & Freeman il989i) . Vallenari et al. ( |2006l l 
found thick disk properties (see their table 6) quite similar to 
the ones obtained in this paper They obtain: crjy=38+7kms"', 
o"f//ci/=1.48, Viag =42 + 7 kms"', and for the thick disk scale 
height: 900 pc. However, they find o-u/o-w-l-9. They also claim 
that "no significant velocity gradient is found in the thick disk", 
implying that the thick disk must be an isothermal component. 



4.5.1. Radial velocities 

The number of RAVE and ELODIE stars used in this analysis is 
a tiny fraction of the total number of stars used from 2MASS 
or UCAC2 catalogues. However they play a key role in con- 
straining Galactic model parameters: the magnitude coverage of 
RAVE stars towards the SGP, from - 8.5 to 11.5, can be 
used to discriminate between the respective contributions from 
each type of star, dwarfs, sub-giants, giants. A future RAVE data 
release (Zwitter et al. submitted) will include gravities, allowing 
for easier identification of dwarfs and red clump giants; it will 
also include element abundances allowing for better description 
of stellar disk populations and new insights into the process of 
their formation. 



5. Conclusion 

We revisit the thin-thick disk transition using star counts and 
kinematic data towards the Galactic poles. Our Galactic model- 
ing of star count, proper motion and radial velocity allows us to 
recover the LF, their kinematic distribution function, their ver- 
tical density distribution, the relative distribution of giants, sub- 
giants and dwarfs, the relative contribution from thin and thick 
disk components, the asymmetric drift coefficient and the solar 
velocity relative to the LSR. 

The double exponential fitting of the vertical disk stellar den- 
sity distribution is not sufficient to fully characterize the thin and 
thick disks. A more complete description of the stellar disk is 
given by its kinematical decomposition. 

From the star counts, we see a sharp transition between the 
thick and thin components. Combining star counts with kine- 
matic data, and applying a model with 20 kinematic components, 
we discover a gap between the vertical velocity dispersions of 
thin disk components with cr^ less than 21 kms"' and a dom- 
inant thick disk component at criy=45.5kms"'. The thick disk 
scale height is found to be 1048+36pc. We identify this thick 
disk with the intermediate metallicity ([Fe/H] — 0.6 to -0.25) 
thick disk described, for instance, by Soubiran et al. (2003). This 
thick disk is also similar to the thick disk measured by Vallenari 
et al ( 120061) who find "no significant velocity gradient" for this 
stellar component. We note that star counts at nij^ ~ 15 suggest 
a second thick disk or halo component with cr^ ~65kms"'. 

Due to the separation of the thin and thick components, 
clearly identified with stars counts and visible within the kine- 
matics, the thick disk measured in this paper cannot be the re- 
sult of dynamical heating of the thin disk by massive molec- 
ular clouds or by spiral arms. We would expect otherwise a 
continuous kinematic distribution function with significant kine- 
matic components covering without discontinuity the range of 
cTw from 10 to 45 km s"' . 

We find that, at the solar position, the surface mass density 
of the thick disk is 27% of the surface mass density of the thin 
disk. The thick disk has velocity dispersions cru - 50kms"', 
cr iv=45. 5 km s"', and asymmetric drift V\ag = 33 ± 2kms"'. 
Although clearly separated from the thin disk, this thick compo- 
nent remains a relatively 'cold' thick disk and has characteristics 
that are close to the thin disk properties. This 'cold' and rapidly 
rotating thick disk is similar to the component identified by many 
kinematic studies of the thick disk (see Chiba & Beers 2000 for 
a summary). Its kinematics appear to be different from the thick 
disk stars studied at intermediate latitudes in pencil beam sur- 
veys (eg Gilmore et al '2002), which appear to be significantly 
affected by a substantial stellar stream with a large lag velocity. 
They interpret this stellar stream as the possible debris of an ac- 
creted satellite (Gilmore 120021 Wyse et al. 2006) . Maybe some 
connections exist with streams identified in the solar neighbor- 
hood as the Arcturus stream (Navarro et al 20041) . 

Some mechanisms of formation connecting a thin and a thick 
components are compatible with our findings. It may be, for in- 
stance a 'puffed-up' thick disk, i.e. an earlier thin disk puffed up 
by the accretion of a satellite (Quinn et al. 1993 ). Another pos- 
sibility, within the monolithic collapse scenario, is a thick disk 
formed from gas with a large vertical scale height before the final 
collapse of the gas in a thin disk, i.e. a 'created on the spot' thick 
disk. We also notice the Samland (2004) scenario: a chemody- 
namical model of formation of a disk galaxy within a growing 
dark halo that provides both a 'cold' thick disk and a metal-poor 
'hot' thick disk. 
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A popular scenario is the 'accreted' thick disk formed from 
the accretion of satellites. If the thick disk results from the ac- 
cretion of just a single satellite, with a fifth of the mass of the 
Galactic disk, this has been certainly a major event in the history 
of the Galaxy, and it is hard to believe that the thin disk could 
have survived this upheaval. 

Finally, from the thick disk properties identified in this paper, 
we can reject the most improbable scenario of formation: the one 
of type 'heated' thick disk (by molecular clouds or spiral arms). 
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